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We calculate the thermal Casimir-Polder potential of Ceo and C70 fullerene molecules near gold 
and silicon nitride surfaces, motivated by their relevance for molecular matter wave interference 
experiments. We obtain the coefficients governing the asymptotic power laws of the interaction in 
the thermal, retarded and nonretarded distance regimes and evaluate the full potential numerically. 
The interaction is found to be dominated by electronic transitions, and hence independent of the 
internal temperature of the molecules. The contributions from phonon transitions, which are affected 
by the molecular temperature, give rise to only a small correction. Moreover, we find that the sizeable 
molecular line widths of thermal fullerenes may modify the nonretarded interaction, depending on 
the model used. Detailed measurements of the nonretarded potential of fullerene thus allow one to 
distinguish between different theories of incorporating damping. 

PACS numbers: 31.30.jh, 12.2Q.-m, 34.35. +a, 42.50.Nn 



I. INTRODUCTION 

It is a remarkable feature of fullerene buckyballs that 
they can exist as delocalized quantum waves as 
is proven almost routinely in matter wave interference 
experiments Q. Such interferometer setups involve 
nanomechanical grating structures, typically made of 
gold or silicon nitride. As the molecules pass through 
the grating slits they experience an attractive dispersion 
force between the polarisable molecule and the grating 
wall. Even though this Casimir-Polder (CP) interaction 
is weak, it must be accounted for in predictions of the 
interference fringes p| 0] . 

The influence of dispersion forces is particularly strong 
in modern near-field interference setups where many dif- 
ferent interference orders contribute resonantly, implying 
that even tiny distortions_of the molecular wave fronts 
affect the fringe pattern In these experiments, it is 
the presence of the Casimir-Polder interaction which im- 
pedes the demonstration of interference with even larger 
and more polarisable particles [§]. At the same time, this 
strong sensitivity of the fullerene matter waves provides a 
means of verifying the precise value and functional form 
of the dispersion forces. It is therefore important to have 
a reliable description of the expected Casimir-Polder po- 
tential available, which should also account for molecules 
not in thermal equilibrium with their environment, given 
that the beam is usually produced by thermal sublima- 
tion. 

Casimir-Polder interactions have been studied inten- 
sively in recent years, though mainly focused on atoms 
[f|. Such studies include thermal equilibrium @ as well 
as non-equilibrium situations @, Q in which the inter- 



nal temperature of the microscopic (atomic or molecular) 
system can be vastly different from that of the macro- 
scopic environment. Carbon-based nanostructures have 
been of particular interest due to applications. In this 
context, the interaction of a carbon nanotube with a sur- 
face has been studied [Toj . 

Motivated by the mentioned matter-wave interference 
experiments we now discuss and evaluate in detail the 
CP interaction of Cgo and C70 fullerenes with planar sur- 
faces made of gold or silicon nitride (SiN^ or S13N4). In 
Sect. [Til we summarise the basic equations that govern 
the CP potential, determine the molecular polarisabilities 
of the fullerenes from spectroscopic data, and list the ma- 
terial parameters of the surface materials. In Sect. IIIH we 
calculate the CP potentials by both numerical and ana- 
lytical means and discuss our results. A short summary 
is given in Sect. IIV1 



II. BASIC EQUATIONS AND PARAMETERS 

We begin by presenting the theory of the thermal CP 
potential and by recording the molecular and material 
properties as obtained from optical data. 



A. Thermal Casimir-Polder potential 

We consider a non-magnetic, isotropic molecule of in- 
ternal temperature T m placed at a distance z from a plane 
non-magnetic surface of permittivity e(to), with both sur- 
face and environment being held at uniform temperature 
T. As shown in Refs. @ and 11 1, the thermal CP poten- 
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tial of the molecule can be given as a sum of non-resonant 
and resonant contributions, 



U(z) = U nTes (z) + U Tes (z). 



(2.1) 



The non-resonant contribution is due to virtual pho- 
tons and it is given by a sum 



U nTes (z) = t^ML [a Tm (i£j) + a Tm (-i£j)] 
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x / d^e- 2re ^[^r s fe,^) 

-(2^ 2 C 2 -^)r p fe, K ± )] (2.2) 

over the purely imaginary Matsubara frequencies i£j with 
£j = (2irkBT/H)j , where the prime indicates that the 
j = term carries half-weight. The properties of the 
molecule are represented by its thermal polarisability 



(2.3) 



with 



71t(w) 



3 huj/(k B T) _ i 



(2.8) 



denoting the thermal photon number. The reflection co- 
efficients for real frequencies can be given as 



r s (u,fc") = , 7>(w,fc ") = —TTT— TX 



(2.9) 



with fc-L^^wV^-fcll 2 if fcH<q;/c fc- 1 = i v /fc ||2 -w 2 /c 2 
if fcH >cj/c and fc^ = y/e{ui)w 2 /c 2 — kW 2 with Im fcj 1 >0. 

We stress that the molecule is out of thermal equilib- 
rium with its environment whenever T m ^ T. Here, the 
thermal polarisability as well as the internal-state pop- 
ulations of the molecule depend on its internal temper- 
ature T m> while the Matsubara frequencies and thermal 
photon numbers are given in terms of the environment 
temperature T. 



where 



Pn(T m ) 



-E n /(k B T ID 



-E k /(k B T m 



(2.4) 



denotes the populations of the molecular eigenstates with 
energies E n and 



i(oj) = lim — 



Wkn I d n k 1 2 



2 n -w 2 -i w (r„ + r fc )/2 



(2.5) 



[ujk n = (Ey. — E n )/h: molecular transition frequen- 
cies, dnk'. electric dipole matrix elements, T n : level 
widths/damping constants, r = 0] are the associated 
polarisabilities. The material properties of the surface 
enter via the reflection coefficients for s- and p-polarised 
waves 



e(i^)/t- L + Ki 

(2.6) 



with = ^K ±2 + [e(i£)-l]£ 2 /c 2 . 

The resonant contribution to the potential is due to 
the absorption and stimulated emission of real photons, 
it reads 



U r cs(z) 



^2p n (T m )l y^[n T (^nfe) + 1] - ^ n T (uj kn ) \ 

n Kk<n k>n ) 

^dfcll^jlm^^d^Ull)] 
^-l)lm(e 2ifc " z r p (K n |,fcll)]l (2.7) 
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B. Molecular and material properties 

The dielectric permittivity e(w) of thin fullerene films 
in the optical regime has been measured by means of 
electron energy-loss spectroscopy [l2| as well as in the gas 
phase [l3| . The frequencies and dipole matrix elements of 
the corresponding electronic transitions of single fullerene 
molecules can be deduced from the data in Ref. in a 
two-step procedure. 

We first apply a simultaneous least-squares fit of an 
rt-oscillator model 



17; LO 



(2.10) 



to the measured data for the real and imaginary parts of 
the permittivity as taken from Ref. [l2| . The fits with 
n = 9 oscillators for Ceo and n — 7 for C70 are illustrated 
in Figs. [T]and [2] The obtained fit parameters Oj, fa and 
r, are given in Tabs. HI and HI1 

Next, we relate the fitted permittivity to the molecular 
polarisability by means of the Clausius-Mosotti law 



a{uS) 



3eo e{uj) — 1 
r) e(uj) + 2 



(2.11) 



where 77 is the number density of fullerene molecules in 
the film. For a face-centred-cubic crystal structure of the 
fullerene molecules in the thin film, one has r\ = 4/a 3 
with a lattice constants a = 1.42 x 10 _9 m for Cgo and 
a= 1.51 x 10 _9 m for C70 [12l ]. We employ a decomposition 
of the resulting expression into partial fractions to write 
it in the form 



a (ui) 
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(2.12) 
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FIG. 1: 9-oscillator fit (solid line) to the spectroscopic data 
(dashed line) of the permittivity of a thin Ceo film. 
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FIG. 2: 7-oscillator fit (solid line) to the spectroscopic data 
(dashed line) of the permittivity of a thin C70 film. 



fli [rad/s] 


fi 


7< [rad/s] 


4.10 x 10 15 


0.120 


5.44 x 10 14 


5.47 x 10 15 


0.663 


1.14 x 10 15 


6.99 x 10 15 


0.664 


1.28 x 10 15 


8.51 x 10 15 


0.348 


1.60 x 10 15 


1.35 x 10 16 


0.0270 


1.72 x 10 15 


1.52 x 10 16 


0.0471 


1.33 x 10 15 


1.85 x 10 16 


0.554 


6.29 x 10 15 


2.54 x 10 16 


0.403 


9.28 x 10 15 


3.27 x 10 16 


0.229 


9.31 x 10 15 



TABLE I: Fit parameters obtained for the permittivity of 
C60 films in the optical regime. Eoo = 1.0463 



fi; [rad/s] 


fi 


ji [rad/s] 


3.76 x 10 15 


0.245 


9.91 x 10 14 


4.97 x 10 15 


0.170 


1.70 x 10 15 


7.92 x 10 15 


1.39 


4.05 x 10 15 


1.58 x 10 16 


0.276 


4.05 x 10 15 


1.98 x 10 16 


0.502 


6.17 x 10 15 


2.51 x 10 16 


0.393 


8.19 x 10 15 


3.34 x 10 16 


0.186 


8.45 x 10 15 



TABLE II: Fit parameters obtained for the permittivity of 
C70 films in the optical 1.0827 



frequency regime. The respective permittivity data ob- 
tained from Fourier transform infrared experiments have 
been fitted to a model of the kind given by Eq. (|2 . 10[) 
with the parameters being listed in Tab. [Vj To ob- 
tain the respective molecular polarisability, we make use 
of the Clausius-Mosotti law (|2 . 1 1 1) . As the measurements 
have been performed at room temperature (T m = 300 K) 
where the phonons are excited to a considerable degree, 
we have to employ the thermal polarisability (I2.3p . Ac- 
cording to Eqs. (|2.4[) and (12.5[) . the latter can be written 
in the form 



from which the transition frequencies Wfco, dipole matrix 2 v tuk n \d n k\ 2 

elements dok and excited-state widths Tk can be read off. a T m l,wj — — ^^PnyJ-m) ■ 

They are given in Tabs. Mil and IIVI Note that we have n,k 

used the ground-state polarisability evo(w) rather than _ 2 )tanhf / ^ LL ° kn 

its thermal counterpart This is a good approx- ~ 3^2^ Pn ' c, - nl > an ^ 2fc B T m 
imation for the considered electronic transitions whose 

frequencies are much higher than the thermal frequency x ^fewjjnfej 

knTm/h=3.93x 10 13 rad/s corresponding to the tempera- w\ n - uj 2 - iuj(T n + T k )/2 



(2.13) 



ture of the molecules in the experiment [12j , T m = 300 K 



In addition to the optical transitions, four phonon with p n k{T) =p n (T)+pk(T). Assuming that all observed 
transitions have been identified for Cqq in the infrared phonon transitions are from the ground state (n = 0), we 
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Lj k0 [rad/s] 


dok 


Cm] 


T k [ra 


d/s] 


4.14 x 


10 15 


7.93 


x 


10~ 


30 


1.10 x 


10 15 


5.73 x 


10 15 


2.36 


x 


lo- 


29 


2.32 x 


10 15 


7.43 x 


10 15 


3.58 


X 


ur 


29 


2.65 x 


10 15 


8.97 x 


10 15 


4.93 


X 


io- 


29 


3.23 x 


10 15 


1.36 x 


10 16 


1.09 


X 


io- 


29 


3.47 x 


10 15 


1.53 x 


10 16 


1.65 


X 


10" 


29 


2.73 x 


10 15 


1.98 x 


10 16 


7.28 


X 


io- 


29 


1.28 x 


10 16 


2.70 x 


10 16 


9.42 


X 


io- 


29 


1.83 x 


10 16 


3.43 x 


10 16 


1.11 


X 


10" 


28 


1.85 x 


10 16 



TABLE III: Transition frequencies, dipole matrix elements 
and widths of the electronic transitions of Ceo- 



Wfco [rad/s] 


d k [Cm] 


Y k [rad/s] 


3.83 x 10 15 


1.40 x 10 -29 


2.01 x 10 15 


5.03 x 10 15 


1.50 x 10~ 29 


3.43 x 10 15 


9.04 x 10 15 


6.91 x 10" 29 


8.12 x 10 15 


1.63 x 10 16 


4.03 x 10" 29 


8.28 x 10 15 


2.10 x 10 16 


7.90 x 10" 29 


1.25 x 10 16 


2.69 x 10 16 


1.15 x 10" 28 


1.60 x 10 16 


3.48 x 10 16 


1.13 x 10 -28 


1.68 x 10 16 



TABLE IV: Transition frequencies, dipole matrix elements 
and widths of the electronic transitions of C70. 



have 



^Po/c(Ti n )tanh( 

k 



3fi 



2k B T u 



Wfeo|di 



■OA' 



— ui 2 — iujTk/2 



(2.14) 



J k0 



The transition frequencies and dipole matrix elements 
can then be readily obtained by comparing with the fit 
for Eq. (|2.1ip . with the results being given in Tab. I VII 

We are going to study the interaction of fullcrcnc 
molecules with Au, SiN x and SiaN4. The permittivity 
of Au can be given as 



e(w) = 1 - 



O 2 
"0 
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lJiUJ 



(2.15) 



fii [rad/s] 


fi 


7i [rad/s] 


9.91 x 10 13 


0.024 


4.33 x 10 11 


1.08 x 10 14 


0.007 


6.03 x 10 11 


2.23 x 10 14 


0.0011 


5.46 x 10 11 


2.69 x 10 14 


0.001 


6.40 x 10 11 



TABLE V: Fit parameters for the permittivity of C60 films 
in the infrared regime. 



Wfeo [rad/s] 


d 0k [Cm] 


T fc [rad/s] 


9.95 x 10 13 
1.09 x 10 14 
2.23 x 10 14 
2.69 x 10 14 


1.69 x IO- 30 
1.00 x lO -30 
5.33 x IO- 31 
5.58 x 10" 31 


8.67 x 10 11 
1.21 x 10 12 
1.09 x 10 12 
1.28 x 10 12 



TABLE VI: Transition frequencies, dipole matrix elements 
and widths of the phonon transitions of Ceo- 



measurements. As reported in Ref. [19(, one has 



Ime(w) = Q(w - Q T ) 



[(w 2 - n 2 y 



~{ 2 U) 2 \uJ 



(2.16) 



with parameters T = 3.48 x 10 15 rad/s, Q, = 1.09 x 
10 16 rad/s, / = 1.13xl0 17 rad/s and 7=1.16x 10 16 rad/s. 
The permittivity at imaginary frequencies as required for 
the non-resonant CP potential can be obtained from the 
Kramcrs-Kronig relation 



2 r° 

em = - \ 

7T Jo 



doi 



ulm£(w) 
~J 2 TW 



(2.17) 



In particular, e(0) — 3.87. 

Alternatively, we also consider non-crystalline SisN4. 
The real and imaginary parts of its permittivity have 
been reported over a wide frequency range [20j . We have 
fitted this data with a single-resonance 4-parameter semi- 
quantum model (2lj . 



n 2 



1W7L 



1W7T 



(2.18) 



The fit, as displayed in Fig. [31 yields the parameters — 
2.69 x 10 16 rad/s, tt T = 1-33 x 10 16 rad/s, -y L = 3.05 x 
10 16 rad/s, and 7 T = 6.40 x 10 15 rad/s. This yields a 
static permittivity e(0) = 4.10. 



where the first term is the response of the conduction 
electrons as described by a Drude model [l5[ and the 
Lorentz-type contributions are due to atomic transitions 
(6-oscillator fit [l6| based on data from Refs. [Tt} and 
18]). Values for the model parameters as taken from the 
above mentioned references are listed in Tab. IVII1 

Diffraction experiments frequently use gratings made 
of low-pressure chemical vapour deposited silicon nitride. 
The imaginary part of the permittivity of this (non- 
stoichiometric) SiN^ has been determined from optical 



III. CASIMIR-POLDER POTENTIAL OF 
FULLERENE 

Using the basic formulas from Sect. Ill Al together with 
the molecular and material parameters from Sec. Ill Bl we 
can now evaluate the CP potential. We begin by calcu- 
lating the CP potential associated with electronic transi- 
tions. As seen from Tabs. HlTl and llVI the electronic tran- 
sition frequencies of Cgo and C70 are much larger than the 
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Q,i [rad/s] 


fi 


7i [rad/s] 





1.37 x 10 16 




5.32 x 10 13 


1 


4.63 x 10 15 


0.762 


1.14 x 10 15 


2 


6.30 x 10 15 


2.41 


2.81 x 10 15 


3 


8.20 x 10 15 


0.0926 


1.52 x 10 15 


4 


1.29 x 10 16 


2.14 


1.06 x 10 16 


5 


2.05 x 10 16 


0.244 


9.12 x 10 15 


6 


3.27 x 10 16 


0.670 


1.37 x 10 16 



TABLE VII: Model parameters for the permittivity of Au. 
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FIG. 3: 4-parameter semi-quantum model fit (solid line) to 
the spectroscopic data (dashed line) of the permittivity Si3N4. 



respective thermal frequency ksT/h = 3.93 x 10 13 rad/s 
even at room temperature. As a consequence, the ther- 
mal photon numbers nT(u> n k) are extremely small. The 
molecules have shown interference at internal tempera- 
tures of 2500 K [23|, [24| , and they are stable up to 6000 K 
_22|. Even at the latter temperature, the thermal fre- 
quency fcB?m/fi = 7.86 x 10 14 rad/s is much smaller than 
the molecular transition frequencies, so the molecule is 
essentially in its electronic ground state, p n (T m ) = 5 n Q. 
As a result, the resonant CP potential (I2.7[) vanishes and 
the CP potential (I2.1[) is entirely given by the nonreso- 
nant contribution (12.21) which simplifies to 



d^e" 2 "^ 



1 TV 



(3-1) 



As first shown in Ref. 25J, we may distinguish three 
asymptotic regimes where the potential reduces to sim- 
ple power laws. At distances much larger than the 



wavelength of the predominant thermal photons, z 3> 
Hc/(k&T), the Matsubara sum is dominated by its first 
term. Using Eq. (|2.6j) and carrying out the k -integral, 
we have 



z 6 



with [ao(0) = ao] 



CzT — 



k B Ta a e(0) - 1 
167nr e(0) + 1 



(3.2) 



(3.3) 



For smaller distances zCfic/^BT); the Matsubara sum 
is well approximated by an integral, so that 



U(z) 



16tt 2 

oo 



poo 

/ da 2 [a (iO + «o(-iO] 
Jo 



1 r„ 



(3.4) 



This distance region can be further divided into the re- 
tarded and nonretarded regimes. At retarded distances 
c/ujka «z« hc/{k B T), the approximations arj(i£) — «o 
and e(i£)~e(0) = £ lead to 



(3.5) 



with 



Ca 



3hcao 
647r 2 e J 1 



dv 



ev - 



ev + \Je — 1 + v 2 

(3.6) 



1 v — \Je — 1 + v 2 
v 4 v + \/e - 1 + v 2 



where we have introduced the new integration variable 
v = k c/£. An explicit formula for C4 and some of its 
limits are found in App. [5] The values of ao are 9.72 x 
10- 39 C 2 m 2 /J and 1.19 x 10- 38 C 2 m 2 /J for C 60 and C 70 , 
respectively. For nonretarded distances z <C c/w/co, an 
asymptotic expansion in terms of z leads to 



with 



C* = 



32tt 2 



•0 Jo 



d£ [a (iO+"oHO] 



gOj) - 1 

e(iO + 1 



(3.7) 



(3.8) 



Using the parameters of Sect. Ill B[ we have calculated 
the values of the coefficients C3, C4 and C3T for Ceo 
and C70 molecules interacting with a perfectly conduct- 
ing surface as well as gold and silicon nitride surfaces. 
The results are given in Tab. IVIIII Due to its larger 
dipole moments, all coefficients are larger for C70 than 
they are for Cgo- The difference is most pronounced in 
the nonretarded regime. Comparing the coefficients for 
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Coefficient — > 


C 3 [Jm 3 ] 


Ca [Jm 4 ] 




Jm 3 ] 


Material 4- 


C60 


C70 


Ceo 


C70 


Ceo 


C70 


Perfect conductor 
Au 
Si 3 N 4 
SiN* 


2.4 x 10" 47 
1.0 x 10~ 47 
8.4 x 1(T 48 
6.3 x 1CT 48 


3.0 x 10" 47 
1.3 x 10~ 47 

1.1 x 10~ 47 
7.9 x 10~ 48 


3.3 x 10" 55 

3.3 x 10 -55 
1.5 x 10~ 55 

1.4 x 10~ 55 


4.0 x 10" 55 
4.0 x 10~ 55 
1.9 x 10~ 55 
1.8 x 10~ 55 


9.0 x 10" 50 
9.0 x 10" 50 
5.5 x 10 -50 
5.3 x 10" 50 


1.1 x 10~ 49 
1.1 x IO- 49 
6.7 x 10" 50 
6.5 x 10~ 50 



TABLE VIII: Coefficients for the asymptotic power laws of the CP potential of fullerene (T = 300K). 
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FIG. 4: Dependence of the retarded CP coefficient C4 of C60 
on the static permittivity. 



the different materials, we note that Au can be consid- 
ered a perfect metal in the retarded and thermal regimes, 
but corrections due to finite reflectivity are quite signif- 
icant at nonretarded distances, leading to a reduction 
of C3 by more than a factor of 2. This is due to the 
fact that the molecular transition frequencies are a size- 
able fraction of the Au plasma frequency. The silicon ni- 
tride potentials are smaller than those of Au due to their 
smaller permittivity. This difference is most pronounced 
at large distances. Comparing the two different silicon 
nitride species reveals that the coefficients for SisN4 are 
larger than those of SiN^ by up to 40% in the nonretarded 
regime. 

The dependence of the retarded CP coefficient C4 on 
the static permittivity of the surface material is displayed 
in Fig.|U The figure shows that the close similarity of the 
C4 coefficients for the two silicon nitride species is due to 
their similar static permittivities. It also reveals that the 
asymptotic C4 value of a metal (e— too) is only reached 
for very large permittivities, hence the large difference 
compared to Au. 

The full potential of Ceo in front of an Au surface has 
been calculated numerically and is displayed in Fig. [5] 
The tabulated results can also be found in the supple- 
mentary materials. The figure shows that the potential is 
faithfully represented by the asymptotic power laws (|3.2I) , 
(I3.5[) and (|3.7[) over a large part of the displayed distance 
range (as indicated by the shaded areas). However, there 



is a large gap between the nonretarded and retarded re- 
gions (between 10~ 7 m and 2 x 10 _8 m) where neither 
limit applies. The potentials for C70 and for different 
surfaces show a similar qualitative behaviour. 

The environment temperature T affects the CP po- 
tential at distances larger than the thermal wavelength 
where thermal photons lead to softening of the potential 
decay. This temperature-dependence is demonstrated in 
Fig. [6] where we display the CP potential for different en- 
vironment temperatures. As seen, thermal photons begin 
to affect the potential at smaller distances for higher en- 
vironment temperatures, resulting in larger long-distance 
potentials. The potentials for 300 K and 600 K begin to 
differ from the zero-temperature result at distances larger 
than about 2 /jm or 4 /jm, respectively. 

In Fig. [JJ we compare the potentials of the two different 
types of fullerenes. The C70 potential is larger than that 
of Cgo by a factor which ranges between 1.3 at small dis- 
tances and 1.2 at large distances. This difference is due 
to the larger dipole moments of C70, recall Eqs. (|3.8[) 
and (|3.6[) for the C3 and C4 coefficients and (|2.5I) for the 
atomic polarisability. The ratio of the potentials drops 
in the transition region between the nonretarded and re- 
tarded regimes. This is because the transition frequencies 
of C70 are slightly larger than those of Cgo- 

The potentials of the two silicon nitride species are 
compared to that of Au in Fig. [8] Note that the curves for 
the two different fullerene molecules are indistinguishable 
in this plot. We see that the ratio is roughly 60% for both 
silicon nitride species at large distances. This similarity is 
due to their very similar static permittivities. At smaller 
distances, the potential of SisN4 is larger than that of 
SiN x by 30%, because its reflection coefficient falls off 
less rapidly with frequency. 

As seen from Tables [TTT1 and HVl the line widths of the 
optical transitions are almost comparable to the transi- 
tion frequencies. In our theory, the impact of the line 
widths is accounted for microscopically (2(|: First, one 
explicitly solves the internal dynamics of the molecules, 
which depends on the transition frequencies and line 
widths. In a second step, this solution is used to deter- 
mine the CP potential (|3.ip , where the dependence on the 
molecular parameters can be expressed via the molecular 
polarisability in a symmetrised form, i[ao(i£) +ao(— 
[26| . As seen from Eq. (|2.5p . the line widths affect the 
imaginary-frequency polarisability most strongly at large 
frequencies. The largest impact on the CP potential 
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FIG. 5: CP potential of Ceo in front of an Au surface at room temperature and its asymptotes. 
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FIG. 7: Comparison of the CP potentials of Ceo and C70 in 
front of gold and silicon nitride surfaces at room temperature. 



FIG. 6: CP potential of C60 in front of an Au surface at 
different temperatures. 

is hence expected at small distances. In Table IVIII[ 
we compare the C3 coefficients including the finite line 
widths with those one would obtain for zero line widths. 
We note that the line widths have practically no influence 
on the nonretarded potential, with differences of about 
1%. This tiny change is quadratic in the line widths and 
unobservable in an experiment. 

In an alternative approach based on linear-response 
theory 0, the atomic properties enter via the 
fluctuation-dissipation theorem. The resulting CP po- 
tential depends on the atomic polarisability in its un- 
symmetrised form, 

h [°° e(ifl - 1 

C 3 ,lrt = „ 2 / dgctoOO ,. c , - -1 ■ (3.9) 




10"' io" B 10" 7 10" 6 10" 5 10" 
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FIG. 8: Comparison of the CP potentials of fullerene in front 
of SiaN4 and SiN x surfaces at room temperature, normalised 
with respect to the corresponding Au potentials. 



The atomic transition frequencies and line widths are not 
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Molecules — > 


Ceo 


C70 


Material 4- 


c 3 


£?3,r fe ->o 


C3,LRT 


c 3 


C3,r fc ->o 


C3XRT 


Perfect conductor 
Au 
Si 3 N 4 
SiN* 


2.36 x 10" 47 
1.01 x 10" 47 
8.45 x 10" 48 
6.26 x 10" 48 


2.34 x 10" 47 
1.00 x 10~ 47 
8.36 x 10~ 48 
6.21 x 10" 48 


2.15 x 10" 47 
9.28 x 10" 48 
7.69 x 10" 48 
5.75 x 10" 48 


2.96 x 10" 47 
1.27 x 10" 47 
1.06 x 10~ 47 
7.86 x 10" 48 


2.93 x 10" 47 
1.25 x 10" 47 
1.05 x 10 -47 
7.76 x 10" 48 


2.68 x 10" 47 
1.15 x 10~ 47 
9.55 x 10~ 48 
7.11 x 10" 48 



TABLE IX: Impact of absorption on the C3 coefficients [Jm 3 ] for the CP potential of fullerene. 



considered explicitly, but only appear at the end of the 
calculation when specifying the polarisability. In par- 
ticular, the line widths now affect the potential already 
to linear order. As seen from Table IIX| this leads to a 
prediction of a reduction of the C3 coefficient by about 
10% due to absorption. This effect could be visible in 
sufficiently accurate C3 experiments, making it possible 
to distinguish between the macroscopic, linear response 
model for absorption and our microscopic model (which 
predicts that the effect of absorption on ground-state po- 
tentials is negligible). This possibility makes fullerenes 
most attractive for CP-potential studies. In contrast, 
atomic systems are unable to resolve the difference be- 
tween the effects of symmetrised vs unsymmetrised po- 
larisabilities (or the neglect of line widths altogether). 

Finally, let us discuss the impact of the infrared res- 
onances on the Ceo potential. As seen from Table I VII 
their dipole moments are much smaller than those of the 
optical transitions (Table HTT|) . On the other hand, their 
transition wavelengths being much longer than those of 
the optical transitions, the nonretarded limit applies over 
a larger range of distances. In addition, the thermal pho- 
ton numbers can take large values even at room temper- 
ature, so that resonant potentials (12.71) come into play. 

At distances up to about 10 fxm, the CP potential due 
to phonon resonances is strongly nonretarded. As shown 
in Ref. (27j, the potential in this regime is well approxi- 
mated by 

tt I \ 3i Phonon /„ 1 n \ 

t'phonon(z) = 5 (3.10) 

with 

C3. Phonon = 77: Pn (T m ) | d-nk | 2 (3-H) 

487T£ 

n.k 

where the sum only runs over phonon transitions. This 
result holds regardless of the environment temperature 
for all materials of sufficiently large permittivity. The 
corresponding C3 coefficient depends on the internal tem- 
perature of the molecule, it ranges from C3,phonon = 
3.4 x 10~ 51 Jm 3 at zero temperature to C3,p n0 non = 2.6 x 
10 _51 Jm 3 at T m = 300K. A comparison with the C 3 co- 
efficients listed in Table I Villi reveals that the potential 
contribution from infrared phonon transitions is smaller 
than the discussed potential from optical transitions by 
more than two orders of magnitude. 



At larger distances, corrections due to imperfect re- 
flectivity manifest themselves [28j . However, they do not 
affect the order of magnitude of the phonon CP poten- 
tial, which remains insignificant. As we have numerically 
verified, the phonon contributions only become relevant 
at very large distances, well beyond 100 /ini. 



IV. SUMMARY 

We have determined the Casimir-Polder interaction of 
Ceo &ud C70 with plane surfaces of Au and two different 
nitride species, as commonly used in molecular interfer- 
ence experiments. The numerically calculated potentials 
is well approximated by the thermal, retarded and nonre- 
tarded asymptotes for large, intermediate and small dis- 
tances, respectively. We have found that the potential is 
entirely due to optical transitions and hence independent 
of the internal temperature of the molecules. The envi- 
ronment temperature affects the potential for distances 
larger than 2 /im at room temperature. Comparing the 
potentials for different silicon nitride species, we have 
found differences of up to 30% in the nonretarded regime, 
which is most relevant for diffraction experiments. 

According to our microscopic theory of the molecule- 
field interaction, the finite line widths of the molecules 
have practically no influence on the Casimir-Polder po- 
tential. A macroscopic linear-response approach, on the 
other hand, predicts that they decrease the nonretarded 
potential by about 10%. The relatively large line widths 
of the fullerene transitions thus make them an ideal sys- 
tem to study the impact of molecular absorption on dis- 
persion interactions. 
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Appendix A: Coefficient for the retarded 
Casimir— Polder potential 



For large and small values of e, C4 behaves like 



The integral of the coefficient C4, given in Eq. Q3.6 
can be solved explicitly with the result 



C 4 = 



3hca f 10 - 3y/e -Ae- 3e 3 / 2 + 6e 2 



647T 2 £ 



3(e-l) 



log 



Ve+T- 1 



- + 21og(Ve + Ve + f) 



2e A - 4e 2 + 3e + 1 
(e - I) 3 / 2 



log( v / i+ Ve- 1] 



(Al) 



<?4 
C 4 



3fea / _ 44 

647r 2 e V 2 \/e 15e 
3ftca /23 169 2 

647r 2 e V30 X ~ 420 X 4 



with x — £ ~ 1 • 



e>l,(A2) 
X«l (A3) 
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